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Simulation of Ground- 
Motion Data 
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(Manuscript received July 31, 1968) 

The time variations of the root-mean-square accelerations, the auto- 
correlation junctions, and the -power spectral density junctions oj 12 strong- 
motion earthquake accelerograms are analyzed. The results indicate that: 
(i) strong-motion accelerograms of sufficiently long duration are stationary 
in the rms sense, (ii) the stationary rms acceleration is a good measurement 
oj earthquake intensity, and (Hi) the autocorrelation and power spectral 
density junctions oj strong-motion accelerograms resemble those oj a nar- 
rowband process. Based upon these results, a method oj determining the 
transfer characteristics of a site is introduced. A procedure for generating 
a filtered, gaussian stationary process to simulate ground motions is de- 
veloped, and two applications of this simulation procedure illustrate its 
significance. 

I. INTRODUCTION 

For many years structural engineers have been concerned with the 
dynamic response of structural systems subjected to seismic excita- 
tions. Ground motions may be caused by natural earthquakes, by 
underground explosions, or by nuclear air blasts. Structures such as 
high-rise buildings, nuclear reactor facilities, or sensitive equipment 
in the vicinity of such events are vulnerable to induced random-type 
disturbances. Traditionally, deterministic methods of analysis relying 
on the known earthquake response spectra have been used. 1 These 
methods have provided valuable information regarding the behavior 
of structures during earthquakes. However, this procedure has a 
serious restriction in that only a few strong-motion accelerograms 
exist which provide ground-motion input. An earthquake is usually 
initiated by a series of irregular slippages along faults, followed by 
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innumerable random reflections, refractions, dispersions, and attenua- 
tions of the seismic waves within the complex ground formations 
through which they travel. 

Since ground motions are generally random, a probabilistic method 
of analysis appears to be more appropriate than the traditional 
method of establishing a reliable design basis for structures sub- 
jected to ground motions. The simulation of ground motion is un- 
doubtedly a necessary step in performing such a probabilistic analysis. 
Because earthquakes are unpredictable, some researchers in structural 
and earthquake engineering have attempted in recent years to simu- 
late earthquakes by using stochastic processes. Both stationary and 
nonstationary models have been investigated. 2-8 

It is the purpose of this paper to investigate the best characteriza- 
tion of ground motions and to establish a valid basis for the stochastic 
simulation of seismic records. For this purpose, 12 commonly used 
strong-motion earthquake accelerograms are analyzed. The time varia- 
tions of the rms accelerations, the autocorrelation functions, and the 
power spectral density functions of these accelerograms are inves- 
tigated. The stationary rms accelerations are used as a measure of 
earthquake intensities and are compared with those found by Housner. 9 
From the power spectral density analysis of existing ground-motion 
records, a linear filter can be determined to represent the transfer 
characteristics of the ground layers at a particular site. This filter 
is used in developing a method of generating a gaussian stationary 
process to simulate ground-motion accelerations. 

In Section III the generation of synthetic ground acceleration 
records using a digital computer is discussed. The synthetic records 
are generated from existing records and from estimated response 
spectra. Two practical examples, of importance to structural engineer- 
ing, are illustrated in Section IV. 

II. STATISTICAL ANALYSIS OF GROUND-MOTION DATA 

In general, no individual record is representative of any other record 
in an ensemble. However, if the data are stationary, valuable statistics 
may be derived by averaging the existing records. If an ergodic process 
is considered, a single record will be sufficient to represent the entire 

process. 

The mean value of a given time-history record x(t) of duration 

T is defined by 

{x)„ = \f *{t)dt. (1) 
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The mean square value of x{t) is defined by 

<*■>« = \ JT* ad dt. (2) 

Following this definition, the root-mean-square or the rms value 
of x(t) is the positive square root of the mean square value (x 2 ) av . This 
is given by 

of 40 = [j, J* x\t) dt^ (3) 

If x(t) is a stationary random process with zero mean, its autocor- 
relation function R x (t) and power spectral density function S x (a>) 
are given by the following transform pair: 

R,(t) = lim | f x(t)x(t + r) dt. (4) 

T— oo J- •'0 

&(«) = j- r R x (r)e- iuT dr. (5) 

Notice that #.,.(t) is always a real-value even function with a 
maximum at t = and that, if (x) av = 0, 

R x (0) = f sju) du> = lim (x\ v , (6) 

•'-oo r-oo 

that is, the maximum autocorrelation represents the mean square 
value of a stationary random process. 

These simple statistical quantities defined in equations 1 through 
5 are useful in applying random vibration theory in earthquake 
engineering. The mean value, mean square value, and root-mean- 
square value represent the time-average strength of the input func- 
tion. The timt variations of (x) av , (x 2 ) nv , or rms of x(t) can be used 
to test the restricted sense stationarity of a time series. The auto- 
correlation function and the power spectral density function are closely 
related to the second-order properties and are generally used as char- 
acterization functions of a stationary process. Since strong-motion 
earthquake accelerograms are, in general, gaussian, 5 either the auto- 
correlation functions or the power spectral density functions will be 
sufficient to provide a complete statistical description. These two func- 
tions also provide a mathematical basis for the random response anal- 
ysis of linear structural systems. 

If x(t) is given in digitized form, equations 1 through 5 can be 
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written respectively as follows: 

1 v 

<Z>av = jy Z) Xk , (7) 

rms of x(t) = [^ E x? J , (9) 



AT-t 



/2 t = R x (k At) = — ^r i; x t x,» , A = 0, 1, 2, • • • , m, (10) 

and 

£ fc = &( w ) =-k + 2E 72,- cos =* + (-1)X,1 , 

fc = 0, 1,2, ■•• , m, (11) 

in which iV = T/At equals the digitization time interval At, x u = 
x(kAt), and m represents the maximum lag number. 

2.1 5ms Strong-Motion Earthquakes 

Using equations 7 through 11, six commonly used strong-motion 
earthquakes each with two horizontal components are analyzed. They 
are: 

A. El Centro, California, December 30, 1934, N-S. 

B. The same, but E-W. 

C. El Centro, California, May 18, 1940, N-S. 

D. The same, but E-W. 

E. Olympia, "Washington, April 13, 1949, S10E. 

F. The same, but S80W. 

G. Taft, California, July 21, 1952, S21W. 
H. The same, but N69W. 

I. Golden Gate Park, San Francisco, March 22, 1957, N10E. 
J. The same, but S80E. 

K. Alameda Park, Mexico City, May 11, 1962, N10°46'W. 
L. The same, but N79°14 , E. 

The accelerograms for these earthquakes are presented in Fig. 1. 
Upon first inspecting these existing strong-motion accelerograms, one 
might conclude that they are nonstationary. However, this conclusion 
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is debatable, considering the lack of sufficient ground-motion data 
for statistical studies, the difficulties of establishing valid nonsta- 
tionary characteristics, and the ultimate objectives in developing a 
useful stochastic model for earthquake-induced ground accelerations. 

The time variations of the rms amplitude of these strong-motion 
accelerograms are obtained and presented in Fig. 2, from which it is 
seen that the earthquake rms acceleration in general approaches a 
stationary value as the duration of the accelerogram is increased. This 
phenomenon indicates that strong-motion accelerograms are sta- 
tionary in their rms amplitude at long duration (T > 20 seconds, ap- 
proximately) . 

The longest and shortest periods existing in a digitized time-history 
record are 2T and 2A£, respectively. These two extremes constitute an 
effective period range for the record. Any analysis of the record be- 
yond its effective period range will be of no significance. A At of 0.01 
second was used in this study. 

The stationary rms acceleration can be used to measure the inten- 
sity of an earthquake if it is accompanied by the corresponding effec- 
tive period range. The rms intensities for all earthquakes, assuming 
T = 20 seconds (corresponding to a period range of to 40 seconds) , 
are listed in Table I. This table shows that the N-S component of 
the El Centro 1940 earthquake (C) has the strongest accelerogram. 

It is interesting to compare earthquake rms intensities with the 
traditionally used Housner intensity, 9 which is defined as 

S/ x = f ° S.(\, T.) dT , (12) 

•'0.1 

where S v is the solution of the following equations: 

u + ^\u + %u= -x 3 {t) (13) 

£,(X, T ) = max | u | (14) 

for constant coefficients A and T . 

Physically, equation 13 represents the equation of motion of a basic 
single-degree-of-freedom linear system with mass m, viscous damping 
A, and stiffness k subjected to earthquake excitation x g {t) at the sup- 
port, as shown in Fig. 3. Correspondingly, T = 2tt {m/k)* is the 
natural period of the system, and S v in equation 14 is the maximum 
relative velocity reached by this system during the earthquake excita- 
tion. Therefore, Housner's intensity definition actually represents an 
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earthquake's potential peak velocity of structures in the period range 
0.1 to 2.5 seconds. In Table II the rras (normalized to a factor of 
2.7/2.20) and Housner's intensities are compared and good agree- 
ments are observed. The basic difference between these two intensity 
definitions is that the rms intensity is independent of the transfer 
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20 30 40 

TIME IN SECONDS 



Fig. 2 — Root-mean-square acceleration of strong-motion earthquakes, (a) 
El Centro and Golden Gate, (b) Olympia and Taft, (c) Alameda Park. 
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U(t) 




X g (t)= GROUND ACCELERATION 



Fin- 3 — Linear mechanical system subjected to ground excitation. 

characteristics of the structure while Housner's intensity depends on 
both the input and the transfer characteristics of the structure. It is 
much simpler to calculate the rms value for a given input time- history 
function than to find the SI\ values, which require numerous mathe- 
matical integrations on the computer. However, the rms intensity 
should be used with care; one should consider the associated effective 
period range and the predominant period contained in the waveform. 

2.2 Determination of the Predominant Period 

The autocorrelation and power spectral density for all accelerograms 
considered are found on a digital computer by using a Fortran program 

Table I — RMS Intensities of Strong-Motion Earthquakes 



Case 


Identification 


RMS acceleration 
at T = 20 b 
(ft/s») 


Average 
rms 

(ft/8 1 ) 


A 
B 

C 
D 


EC34NS 
EC34EW 
EC40NS 
EC40EW 


1.4 

1.4 
2.4 
2.0 


1.4 
2.20 


E 
F 

G 
H 


OL49S10E 
OL49S80W 
TA52S21W 
TA52N69W 


1.7 
2.0 
1.4 

1.4 


1.85 
1.4 


I 

J 
K 
L 


GG57S80E 
GG57N10E 
AL62N10°46'E 
AL62N70°14'E 


0.8* 
0.6* 
0.75 
0.7 


0.7* 
0.725 



♦Taken at T = 10.0s. 
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based on equations 10 and 11. Those functions, R(t) and S(o)), for the 
El Centro 1940 earthquake are shown in Figs. 4 and 5 for the N-S 
and E-W components, respectively. The smooth curves in the power 
spectral density diagrams (Figs. 4b and 5b) were obtained by introducing 
a Hanning or smoothing procedure 10 to the raw estimates given by 
equation 11. The autocorrelation functions (Figs. 4a and 5a) all have 
a maximum at r = and diminish rapidly at large correlation time. 
The ordinates of the power spectral density functions (Figs. 4b and 5b) 
generally increase with increasing frequencies to a maximum value at 
some frequency which may be considered a predominant or character- 
istic ground frequency, and then decrease rather rapidly toward zero 
in an asymptotic manner. Also of interest is the fact that this general 
rise and fall of the power spectral density function is accompanied by 
local random fluctuations. 

The above results indicate that the autocorrelation and power spec- 
tral density of strong-motion earthquake accelerograms resemble those 
of a narrowband process. This implies that earthquake acceleration 
may be simulated by passing a wideband process through a linear 
filter which reflects the local geological conditions. 

III. STOCHASTIC SIMULATION OF GROUND-MOTION ACCELEROGRAMS 

3.1 Basic Requirements for Ground-Motion Model 

There are many occasions when a ground-motion model is required. 
Examples are the prediction of ground motion at a certain site where 
no past records are available, and statistical analyses of structural 
responses based upon very limited actual ground-motion records. The 
hypothesized models must, of course, possess the pertinent characteris- 
tics of real ground motions and must be supported by existing data. 
More importantly, these models must properly reflect the damage (or 
response) potential of future ground motions to a wide range of struc- 
tures. 



Table II — Comparison of Earthquake Intensities 



Earthquake 
case 


A & B 


C & D 


E&F 


G & H 


IZ& J 


K& L 


Normalized 
rms intensity 


1.7 


2.7 


2.2 


1.7 


0.84 


0.9 


Housner 
intensity 


1.9 


2.7 


1.9 


1.6 


unav 


ailable 
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Fig. 4_N-S component of earthquake at El Centro, Calif., May 18, 1940. 
(a) Autocorrelation function, (b) Power spectral density function, (c) Compari- 
son of actual and idealized power spectral density functions. 
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Fig. 5 — E-W component of earthquake, El Centro, Calif., May 18, HMO. 
(a) Autocorrelation function, (b) Power spectral density function, (c) Compari- 
son of actual and idealized power spectral density functions. 
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The random phenomenon observed in earthquake ground motions 
suggests that one can reasonably assume that the resultant seismic 
wave arriving at the surface of the ground will contain a random as- 
semblage of velocity impulses having gaussian distribution. Further- 
more, the analyses of the earthquake accelerograms in the previous 
sections show that they have characteristics resembling those of a 
narrowband process. It is therefore postulated that the use of stochas- 
tic processes would be appropriate in modeling earthquake ground 
motions. 

Basic criteria for establishing such a stochastic model can be spe- 
cifically stated as follows: 

(i) The model must have the basic properties reflected by the past 
recorded data such as the intensity, duration, general physical ap- 
pearance, and all important characteristics resulting from local geo- 
logical conditions. 

(it) The response statistics of the stochastic model must be equi- 
valent to those produced by the real ground motion or to those 
predicated on strong theoretical or empirical bases. 

Some simplifying assumptions required in the development of our 
stochastic model for earthquakes are: 

(i) The input seismic wave transmited at bedrock by an earth- 
quake is represented by stationary, white noise. 

(U) The ground layers of a seismic station during the shock are 
represented by a single-degree-of-freedom system with linear be- 
havior. (A more sophisticated ground-layer filter may also be used.) 

(Hi) Local random fluctuations appearing in the power spectral den- 
sity function of the real earthquake accelerogram are neglected when 
modeling the transfer characteristics of the site. 

The representation of ground layers by a single-degree-of-freedom 
system is shown in Fig. 6. Such a simple system is characterized by its 
transfer function h(r) in the time domain or H{iw) in the frequency 
domain, with h(r) and H(iu) known as the unit impulse response and 
the complex frequency response of a linear system, respectively. 

If the simple mechanical system shown in Fig. 3 is used as the linear 
filter, the power spectral density of the total acceleration of the mass, 
when the acceleration at the support is taken as the input, can be 
found in terms of the corresponding transfer function. 

S» = I HM [ S< - [x _ {uM y + 4X ; (w/ttp) » W 
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GROUND 
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SEISMIC STATION 
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.BEDROCK 






SHOT (SEISMIC WAVES) 



OUTPUT = GROUND 
ACCELERATION 



|° LINEAR 
SYSTEM 



h(T), H(Loj) 



INPUT = WHITE NOISE 



Fig. 6 — Ground layers represented by a linear system. 

where S { is the constant power spectral density of the input acceleration 
by assumption (i) ; X„ and u B are the ground characteristic damping and 
frequency, respectively. 

For known *$?,• , X„ , and w B , equation 15 gives a smooth curve for 
S (u>). The location of the peak amplitude of S (u) is determined by 
o)„ . The shape or the rate of the rise and fall of S (a)) is governed by 
X„ , and the relative amplitude of S (a>) depends on S { . For any power 
spectral density derived from a real ground-motion record, an idealized, 
smooth equivalent power spectral density can be found by using equa- 
tion 15 and by properly adjusting the three-characteristic constant 
X„ , cd„ , and Si in this equation. The idealized power spectral density 
function for both horizontal components of the El Centro, California 
1940 earthquake were found by the prescribed procedure. As shown in 
Figs. 4c and 5c for both cases, the idealized power spectral density 
function covers the same area as the actual power spectral curve in the 
frequency range through 10 cps. The «, , X, , and £,• values used for 
the N-S and E-W components of the El Centro 1940 earthquake are 
15.5 rad/s, 0.42, 0.046 ft 2 /s 4 /rps, and 14.7 rads/s, 0.41, 0.033 ft 2 /sVrps, 
respectively. 

Once the representative ground-layer filter for a given ground motion 
is determined, the artificial earthquake can be generated by passing 
white noise through this filter and measuring the output time history. 



3.2 Generation of a Gaussian Stationary Process 

This procedure starts by sampling a sequence of pairs of statistically 
independent random numbers X\ , Xi ; x 3 , .t 4 ; • • • ; .r„_i , x n all of which 
have a uniform probability distribution over the range < x < 1. A 
new sequence of pairs of statistically independent random numbers 
Vi i 2/2 ; 1/3 , Vi ; • " • ; y n -\ , y„ arc then generated using the relations 
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y t = (-2 log. a,) 1 cos 27rx i+1 j « i f 3, ...,»- 1 (16) 

y>+i = (-2 log. Xi) h sin27rx, +1 
which have been shown to have a gaussian distribution with a mean 
of zero and a variance of unity. 11-13 

A single waveform y(t) can now be established by assigning the 
values 2/1 , #2 . • • , T/n to n successive ordinates spaced at equal intervals 
At along a time abscissa and by assuming a linear variation of ordi- 
nates over each interval. To define a time origin, assume that the 
initial ordinate y , which is taken equal to zero, is located at t = t 
where t is a random variable having a uniform probability density 
function of intensity 1/At over the interval < t < At. 

For practical reasons At must, of course, be taken as finite; how- 
ever, its value should be set sufficiently small so that the true power 
spectral density function is reasonably constant at intensity S f over 
the lower range of frequencies which are to be properly represented in 
the process. A value of 0.025 second or smaller is recommended. 

To establish the desired stationary process a(t), each member y r (t) 
(r = 1, 2, • • • , N) of the normalized process y{t) must be filtered in 
accordance with equation 15. This step can be accomplished by assum- 
ing that a simple single-degree-of-freedom system having an undamped 
circular frequency o>, and a damping ratio X B is subjected separately to 
the N support accelerations y r (t) and then by calculating the correspond- 
ing absolute or total acceleration functions a T (t) of the mass. In mathe- 
matical form this statement is equivalent to saying that one must solve 
the differential equations 

Z r (t) + 2o>„X,Z r (fl + ufZr® = -y r (t) r = 1, 2, • • • , N (17) 

for the functions 2 r (t) and then evaluate the desired family of accelera- 
tion functions a r (t) using the relation 

o r («) = Z T (t) + Vr(t) r = 1, 2, • • • , N. (18) 

Equation 17 can be solved numerically on a digital computer using 
the standard linear acceleration method. Based on the assumption that 
the input earthquake duration is divided into the very short equal 
time intervals At, and the response acceleration varies linearly over 
each interval, this method gives the following simple relations be- 
tween the response displacement Z r and its derivatives at step n and 
n+1: 
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(ZX At , (Z r ) n+1 



(Z r ) n+1 = (Z,). + ^"- At + 



A/ 



(Z r ) B+1 = (ZX + (Z r )„ A/ + ^s A/ 2 + ^^ A/ 2 . (19) 

•5 fa 

IV. ILLUSTRATIONS OF GENERATION OF ARTIFICIAL ACCELERATION 

4.1 By Known Power Spectral Density 

As the first example, artificial earthquakes are generated using the 
prescribed approach, to simulate the average of U. S. strong-motion 
earthquakes observed at seismic stations having a firm soil founda- 
tion. Values of 15.6 rad/s for w g and 0.6 for k g , representative of such 
soil conditions, were used. A total of 50 artificial accelerograms (N 
= 50) were generated for process a(t) with a duration of 30 seconds 
which corresponds to M = 0.025 second and n = 1200. The intensity- 
Si of the unfiltered "white noise" was set at 0.00814 ft 2 /s 3 so that the 
mean velocity response spectrum curves for the filtered process a(t) 
would give a "best fit" with the standard response spectrum curves 
published by G. Housner. 9 This intensity is slightly less than the 
value of 0.0063 ft 2 /s 3 used by J. Penzien" in a previous investigation 
to correlate the mean velocity response spectrum curves for "white 
noise" with Housner's standard curves. 

A typical sample member of artificial accelerograms generated is 
shown in Fig. 7. It is interesting that this accelerogram appears to be 
very similar to the real accelerograms in Fig. 1 except for the general 




10 15 20 

TIME IN SECONDS 



Fig. 7 — Sample member of artificial accelerogram, 
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0.6 



0.4 - 



K 0.2 H 



-0.2 




1.0 1.5 2.0 

TIME IN SECONDS 



3.0 



Fig. 8 — Autocorrelation function of the sample member of artificial accelero- 
gram. 

stationary apparance of the artificial accelerograms versus the non- 
stationary appearance of the real accelerograms. 

The autocorrelation and corresponding power spectral density func- 
tion for this sample artificial earthquake are shown in Figs. 8 and 9, 
respectively. It is of particular significance that, while these autocor- 
relations and power spectral densities are similar to those for the real 
earthquake (Figs. 4 and 5), local random fluctuations also appear in 



0.016 




10 20 30 40 50 

FREQUENCY IN RAD/S 

Fig. 9 — Power spectral density function of the sample member of artificial 
accelerogram. 
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Fig. 10 — Comparison of velocity spectra of earthquakes, X = 0.02. 

the individual power spectral density function. It is expected that 
these local fluctuations and the variations from one individual power 
spectral density function to another will be eliminated by the averag- 
ing procedure. It has been shown that the average power spectral den- 
sity function of artificial earthquakes is quite close to the prescribed 
function. 15 

The average response spectra (A = 2%) of 50 artificial earthquakes 
are compared with the real spectra given by Housner in Fig. 10. Good 
agreements are observed over the significant period range. The con- 
fidence limits for these average spectra are given in Fig. 11 by extend- 
ing 3cr (standard deviation) of S v above and below the mean. Since 
the distribution of the response spectra is no longer gaussian, these 
limits correspond to an approximate 89 percent confidence band 
according to Chebyshev inequality. "■ 17 

The above results (Figs. 7 through 11) show that the basic require- 
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Fig. 11 — Variations and confidence limits for average velocity spectra of 
artificial earthquakes, X = 0.02. 
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ments for ground-motion simulation have been satisfied. It takes ap- 
proximately 10 seconds to generate one sample function on a CDC 
6600 digital computer. These facts indicate that a good, reasonably 
economical ground-motion simulation has been achieved. 

4.2 By Knovm Response Spectrum 

In the second example, we have a situation somewhat different from 
the first. Engineers are asked to generate representative ground-mo- 
tion accelerations from an estimated pseudovelocity spectra as given 
in Fig. 12. There is no past recorded ground-motion data available. 
The generated artificial earthquakes will be used to evaluate the prob- 
able damage of structures resulting from possible seismic events tak- 




CONFIDENCE _ 

BAND 

CONFIDENCE 



MEAN 



1 2 

PERIOD IN SECONDS 



Fig. 12 — Pseudovelocity response spectra, \ = 0.02. 
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ing place in the vicinity of a proposed site. The estimated response 
spectra were obtained from a careful examination of earthquake ac- 
tivities in a certain seismic area and from an analogous prediction 
using records taken at stations having geological conditions similar or 
comparable to those of the given site. If the event refers to an under- 
ground explosion, the estimated response spectra might be predicted 
from past events under the same physical conditions (yield level and 
depth of the explosion, epicenter distance, and so on) and similar 
geological environments. 

The given pseudoresponse spectrum curves (Fig. 12) show peaks at 
period T. ^ 1.25 seconds if the long-period portion (T > 2.5 seconds) 
or the low-frequency portion (u n < 2.5 rad/s) is neglected. The cor- 
responding frequency value where the peak spectrum amplitudes occur 
is approximately 5 rad/s. This value can be used as an approximate 
predominant ground frequency u, . The proposed site has a firm soil 
foundation, therefore a value of 0.6 for the characteristic ground damp- 
ing X„ is used. 

Using these two characteristic values and letting S t be unity in 
equation 15, the same procedure as used in the first example can be 
followed again. Five typical artificial accelerograms of 30-second 
duration are generated and shown in Fig. 13. The individual response 
spectra and the mean response spectra are obtained (Fig. 14). The 
general shape of the pseudoresponse spectra using Si = 1.0 is similar to 
the estimated spectra. By matching the area covered by the estimated 
and pseudospectra curves, the amplitudes of the pseudospectra are 
normalized to have the same order of magnitude as the estimated 
spectra. The accelerations shown in Fig. 13 which have been modified 
by these same normalization procedures, represent the final form of 
the desired artificial accelerograms. 

The comparisons of individual and mean pseudoresponse spectra 
with the given response spectra (Fig. 14) are satisfactory except for 
the high-period range (T > 2.5 seconds). If only those general building 
structures are considered which have a natural period range of 0.3 to 
2.5 seconds, the large discrepancy of response spectra in the high- 
period range is not significant. The mean pseudoresponse spectrum 
(Fig. 14f) gives better results, as is expected. Overall results of the 
simulation may be improved with more accurate determination of values 
for a, , X„ , and the larger sample size N of the simulating process. 
Methods for such improvements are presented in a separate paper. 18 
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Fig. 13 — Pseudoearthquake accelerograms. 
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V. CONCLUSIONS 

From the results presented in this investigation, the following con- 
clusions may be deduced: 

(i) In general, strong-motion earthquake accelerograms of suf- 
ficiently long duration (T > 20 seconds, approximately) are station- 
ary in the rms sense. 

(u) The rms acceleration of an earthquake, along with the effec- 
tive period range and the predominant period of the associated ac- 
celerogram, can be used to determine the intensity of the earthquake. 

(Hi) The autocorrelation and the power spectral density of a strong- 
motion earthquake resemble those of a narrowband process. 

(iv) The power spectral density analysis of existing earthquake 
records suggests that the transfer characteristics of the given site can 
be represented by a simple, linear system. 

(v) A filtered, gaussian, stationary process generated on a digital 
computer proves to be successful in modeling ground motions in- 
duced by earthquakes or by nuclear underground explosions. 
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